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Abstract. The apparent discrepancy between low and 
high D abundances derived from QSO spectra may be 
caused by spatial correlations in the stochastic velocity 
field. If one accounts for such correlations, one finds good 
agreement between different observations and the theo- 
retical predictions for standard big bang nucleosynthesis 
(SBBN). In particular, we show that the H+D Lya profile 
observed at z a = 0.7 toward Q 1718+4807 is compatible 
with 4.1 10~ 5 < D/H < 4.7 10~ 5 . This result is consistent 
with our previous D/H determination for the z a = 2.504 
system toward Q 1009+2956 and, thus, supports SBBN. 
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1. Introduction 

From recent HST observations of a low-redshift (z a = 0.7) 
absorption-line system toward the quasar Q 1718+4807 
(z e = 1.084) Webb et al. (1997a,b) deduced D/H = 
1.8 - 3.1 10~ 4 . This ratio is significantly higher than that 
derived from other quasar spectra at z a = 2.504 [D/H = 
1.8 - 3.5 10~ 5 by Tytler & Buries, 1996; D/H = 2.9 - 
4.6 10~ 5 by Levshakov, Kegel & Takahara, 1997 (LKT, 
hereinafter)], and at z a = 3.572 [D/H = 1.7 - 2.9 10~ 5 by 
Tytler et al, 1996; D/H > 4 10~ 5 by Songaila et al, 1997]. 

The apparent spread of the D /H values leads some au- 
thors to assume fluctuations in the baryon-to-photon ratio 
at the epoch of BBN (see e.g. Webb et al, and references 
cited therein). On the other hand, according to the ba- 
sic idea of homogeneity and isotropy of big bang theory 
the primordial deuterium abundance should not vary in 
space. One can only expect that the D/H ratio decreases 
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with cosmic time due to conversion of D into 3 He and 
heavier elements in stars. To check whether big bang nu- 
cleosynthesis has occurred homogeneously or not, precise 
measurements of absolute values of D/H at high redshift 
are extremely important. The fundamental character of 
this cosmological test requires an unambiguous interpre- 
tation of spectral observations. 

It is well known, however, that the physical parame- 
ters derived from spectral data depend on the assumptions 
made with respect to the line broadening mechanism. For 
intcrgalactic absorption lines a 'non-thermal broadening' 
is usually assumed to be caused by large scale motions 
of the absorbing gas. The commonly used microturbulent 
approach disregards all correlations of the velocity field, 
implying a symmetrical (Gaussian) distribution of the ve- 
locity components parallel to the line of sight and a sym- 
metrical line profile. 

Actually, any turbulent flow exhibits an immanent 
structure in which the velocities in neighboring volume el- 
ements are correlated with each other. Different aspects of 
the line formation processes in correlated turbulent media 
have been discussed recently in a series of papers by Lev- 
shakov & Kegel (1997, LK hereinafter), Levshakov, Kegel 
& Mazets (1997, LKM hereinafter), and by LKT. 

Once the statistical properties of a turbulent cloud 
have been specified, the way spectral lines ought to be cal- 
culated, depends on the problem considered. Considering 
emission lines one is dealing with many lines of sight and 
the observed intensity should closely correspond to the 
theoretical expectation value (see e.g. Albrecht & Kegel 
1987). If, however, one observes the absorption spectrum 
in the light of a point-like background source, the actually 
observed line profile is determined by the velocity distri- 
bution along the particular line of sight. Therefore, the 
intensity may deviate substantially from the expectation 
value, since averaging along one line of sight corresponds 
to averaging over an incomplete sample (for details see 
LK and LKM). For large values of the ratio of the cloud 
thickness L to the correlation length / the distribution 
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Fig. 1. Observations and RMC fits for Q 1718+4807. a 
HST/GHRS data (open circles) and calculated profiles for Hi 
(dashed curve), Dl (dotted curve) and H+D (solid curve). The 
latter correspond to model (e) in Table 1. The spectral reso- 
lution is 0.1 A (FWHM). b Residuals e in units of a no i S G (see 
text), c IUE spectrum (open circles) and fit (solid curve). The 
spectral resolution is 2.95 A (FWHM) 



function p(v) for the velocity component parallel to the 
line of sight approaches the statistical average, which has 
been assumed to be a Gaussian. For values of L/l of only a 
few, however, p(v) may deviate substantially from a Gaus- 
sian, and is asymmetric in general. This leads to a complex 
shape of the absorption coefficient for which the assump- 
tion of Voigt profiles could be extremely misleading. The 
actual D /H ratio may turn out to be higher or lower than 
the value obtained from the Voigt-fitting procedure. 

The present Letter is primarily aimed at the inverse 
problem in the analysis of the H+D Lya absorption ob- 
served by Webb et al. (1997a, b). The original analysis was 
performed in the framework of the microturbulent model. 
Here we make an attempt to re-analyze the observational 
data on the basis of a more general mesoturbulent model. 
We consider a cloud with a stochastic velocity field with 
finite correlation length but of homogeneous (H I-) density 
and temperature. The velocity field is characterized by its 
rms amplitude a t and its correlation length I. The model 
is identical to that of LKT. - The objective is to investi- 
gate whether the data in question may also be interpreted 
by a lower D/H ratio consistent with the values found for 
other absorption systems. 
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Fig. 2. As Fig. 1 but for model (c) in Table 1 
2. Parameter estimation 

To estimate physical parameters and an appropriate veloc- 
ity field structure along the line of sight, we used a Reverse 
Monte Carlo [RMC] technique (see LKT). The algorithm 
requires to define a simulation box for the 5 physical pa- 
rameters : N(Hi), D/H, Tkin, <Jt/v t h, and L/l (here v t h 
denotes the thermal width of the hydrogen lines). - The 
continuous random function of the coordinate v(s) is rep- 
resented by its sampled values at equal space intervals As, 
i.e. by {vi, v 2 , . . . ,Vk}, the vector of the velocity compo- 
nents parallel to the line of sight at the spatial points Sj . 

In the present study we adopted for the physical pa- 
rameters the following boundaries : 

According to Webb et al., N(Hi) is well restricted within 
the range from 1.70 10 17 cm~ 2 to 1.78 10 17 cm -2 , as de- 
rived from the observed Lyman-limit discontinuity. 

For D/H we use the range from 3.0 10~ 5 to 5.0 10 -5 , trying 
to find a low D/H solution. 

To restrict Tki n and a t /v t h one has to assume a model for 
the absorbing material. It is generally believed that ab- 
sorption line systems with N(Hi) ~ 10 17 cm -2 arise in 
the halos of putative intervening galaxies. Direct observa- 
tions of galactic halos at z > 2 (van Ojik et al. 1997) show 
that a t ~ 40 ± 15 km s" 1 , if T kin ~ 10 4 K. For T kin we use 
the interval 10 4 — 2 10 4 K, and, thus, <7t/«th may range 
within 1.3 - 4.3 . 

For L/I > 1 the meso- and microturbulent profiles tend 
to be identical (see LK). Webb et al. thoroughly investi- 
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Fig. 3. Probability density functions (p.d.f.) for residuals e 
shown by dotted line (corresponds to Fig. lb) and by solid line 
(corresponds to Fig. 2b) and their fit by a Gaussian with mean 
E{e) = and variance Var(e) — 1 



gated different microturbulent models. Therefore we con- 
sider only moderate L/l ratios in the range 1.0 - 5.0. 

Following Webb et al. (1997b), we exclude the Sim 
line from our analysis of the H+D Lya profile. As shown 
by Vidal-Madjar et al. (1996), 'deducting lines of sight ve- 
locity structure for D/H evaluations from ionized species 
could be extremely misleading'. But we fix 2; a (Siin) = 
0.701024 as a more or less arbitrary reference radial ve- 
locity at which Vj =0. 

Having specified the parameter space, we construct the 
objective function (LKT) : 



I(\i)-r(Xi) 



(1) 



where r(Xi) is the simulated random intensity, /(A,) the 
observed normalized intensity within the zth pixel of the 
line profile, cr noisc the experimental error level, and v = 
m — n the degree of freedom (to is the number of data 
points and n is the number of fitted physical parameters, 
n = 5 in our case). 



Here a,. 



is assumed to be constant over the en- 



tire H+D profile. In order to estimate (T no i so we have 
chosen a portion of the Q 1718+4807 continuum from 
the left (AA L = 2066.02 - 2066.90 A) and the right 
(AA R = 2068.52 - 2070.01 A) hand side of the hydrogen 
absorption (see Fig. 1), as well as the bottom of the Lya 
line (AA B = 2067.35 - 2068.30 A). These three regions 
yield ctl = 0.095, or = 0.093, and ob = 0.086, respec- 
tively. At z — 0.7 where the number density of the Lya 
lines is far lower than at high rcdshifts, the experimental 
noise is expected to be free from the contamination by 
weak Lya forest lines. Therefore the estimated noise level 
corresponds to the experimental uncertainties. To approx- 
imate the error level we set a no i 8C — 0.09. 

Finally, our objective function contains only the blue 
and the red wing ofthcHiline (AA = 2066.90-2067.35 A 
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Fig. 4. Probability density functions (p.d.f.) for velocity com- 
ponents p(v) for the RMC solutions shown in Fig. la (dotted 
line histogram) and in Fig. 2a (solid line histogram). Both are 
blue-shifted by about —20 km s _1 with respect to z a (Sim). 
Note the asymmetric shapes. For comparison, the short dashed 
line shows p(v) for the microturbulent model adopted by Webb 
et al. 



and AA = 2068.30-2068.52 A, respectively) since they are 
more sensitive to the parameter variations than the central 
part of the line is. By this we restrict the total number of 
data points to to = 30, implying v ~ 25 and -^xt a — 1-174 
for the credible probability P a = 1 — a = 0.75. 

The estimated parameters for a few adequate RMC 
profile fits (Xmin < X25 0.25) are listed in Table 1. The 
derived deuterium abundance (D/H) ~ 4.4 10~ 5 is about 
4-7 times smaller than the limiting values of 1.8 — 3.1 10 -4 
found by Webb et al. (1997b) in the microturbulent model 
excluding the Si in line. To illustrate our results, we show 
in Figs. 1 and 2 H+D Lya profiles for the two calculations 
with the lowest and the highest D abundances found in the 
mesoturbulent model (D/H = 4.111 10~ 5 and 4.755 10" 5 , 
respectively). They are shown by the solid curve, whereas 
the open circles give the experimental intensities I(Xi). 

The residuals = [/(A;) — r(Ai)]/er n oisc shown in Figs, 
lb and 2b by filled circles are normally distributed with 
zero mean and unit variance (see Fig. 3). This fact, rather 
trivial in case of high S/N data, becomes crucial when 
dealing with spectra as noisy as the present one is, because 
of the probability to be trapped into fitting the noise fea- 
tures. The good concordance of the residuals with the ex- 
pected normal distribution is here a significant argument 
for the validity of the results obtained. 

To check the w(s)-configurations estimated by the 
RMC procedure, we calculated profiles for the higher order 
Lyman lines and the shape of the Lyman- limit discontinu- 
ity and then superposed them to the corresponding part 
of the IUE spectrum. The results are shown in Figs, lc 
and 2c where again open circles correspond to the ob- 
served intensities and the computed spectra are shown by 
the solid curves. We do not find any pronounced discor- 
dance of calculated and real spectra. On the other hand, 
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Table 1. Cloud parameters derived from the Lya profile by the 
RMC method [Ni7(Hl) column density in units of 10 1 cm - , 
D/H in units of 10~ 5 , T^km kinetic temperature in units of 10 4 
K, a t turbulent velocity in km s _1 ] 





Ni 7 (Hi) 


D/H 




at 


L/l 


iy 2 - 

jy A mm 


(a) 


1.732 


4.565 


1.41 


22 


2.7 


1.064 


(b) 


1.739 


4.562 


1.60 


18 


3.9 


1.002 


(c) 


1.759 


4.755 


1.75 


40 


2.8 


1.106 


(d) 


1.761 


4.555 


1.46 


29 


4.3 


1.086 


(e) 


1.768 


4.111 


1.51 


26 


3.5 


0.897 


(f) 


1.771 


4.442 


1.76 


23 


3.4 


1.162 


(g) 


1.776 


4.249 


1.62 


28 


4.0 


1.114 



the spectral resolution of 2.95 A is not sufficient to follow 
a fine velocity field structure within the z a = 0.7 absorber. 

The derived u(s)-configurations are not unique. Ta- 
ble 1 demonstrates the spread of the rms turbulent veloc- 
ities from 18 up to 40 km s -1 . ft is worthwhile to em- 
phasize once more that the projected velocity distribution 
function may differ considerably from a Gaussian. Fig. 4 
shows an example of such distortions caused by a poor sta- 
tistical sample (i.e. incomplete averaging) of the velocity 
field distributions for the two cases of the lowest and the 
highest D/H ratios from Table 1 [model (e) and (c), re- 
spectively]. Both p(v) distributions are asymmetric. This 
is the main reason why the absorption in the blue wing of 
the H I Lya line may be enhanced without any additional 
Hi interloper(s). 

3. Conclusion 

We have shown that the interpretation of the HST and 
IUE spectra obtained by Webb et al. is not unique. The 
data can as well be modeled with a low value of the D /H 
ratio if one accounts for spatial correlations in the large 
scale velocity field. 

The RMC results may be tested, in principle, by addi- 
tional observations of higher order Lyman lines with the 
same spectral resolution as Webb et al. used for Ly-a. 
Indeed, if p(v) is asymmetric, this will show up in the 
profile shapes of the higher order Lyman lines. For the 
physical parameters listed in Table 1, the effect becomes 
visible starting from Ly-4 [the Ly-a, -/?, -7 lines are insen- 
sitive to t,hp asymmetry of p(v) Hue to their high optical 
1 
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Fig. 5. The Hi Ly-4 and Ly-12 mesoturbulent spectra for 
model (c) (solid curves) and (e) (dotted curves) in Table 1. 
Short dashed curves show microturbulent profiles calculated 
for the mean N(Hl) = 1.738 10 17 cm -1 and the mean Doppler 
parameter 6(Hl) = 25.5 km s _1 of the data by Webb et al. The 
spectra are convolved with a Gaussian instrumental profile of 
FWHM = 0.1 A 
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Fig. 4 - dotted curves for model (e) and solid curves for 
model (c) of Table 1. The line shapes clearly depend on 
the velocity field structure and are asymmetric in general. 
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